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Abstract 

A high-school exercise is used to get an insight into planetary motion. 

1 Introduction: some history 

Science contests for high-school students have a hundred-years-old tradition in Hungary ^ |2J . 
The first one, open to students who just completed their high-school cursus, was organized in 
f894 to mark the appointment of L. Eotvos as Minister of Education. This was also the year 
when the High-School Journal for Mathematics and Physics was first published. This Journal 
sets selected problems each month, and the best solutions are published under the name of the 
student who solved it. 

The early prize-winners of the Eotvos contest include such outstanding future scientists as L. 
Fejer (who became famous in mathematical analysis), Th. von Karman (who made important 
contributions to hydro and aerodynamics), A. Haar (remembered for his invariant measure on 
group manifolds) . . . followed by many others. 

In 1916, in the middle of World War I, a physics contest was held for the first time - and 
L. Szilard (who later, with E. Fermi, patented the nuclear reactor) came second. E. Wigner 
(who got the Nobel prize for using group theory in atomic physics), and J. von Neumann (who 
built the first computer and contributed to many fields, ranging from game theory to quantum 
mechanics) only missed the contest because of the turmoil following the disaster of World War I. 

In 1925 Edward Teller not only got the first prize in physics, but also shared the mathematics 
prize with Laszlo Tisza. They became friends and their friendship lasted until Teller's recent 
death. Soon Tisza (then a mathematics student in Gottingen) changed to physics under the 
influence of Born, and published his very first paper jointly with Teller (who was preparing his 
Ph. D. in Leipzig with Heisenberg) on molecular spectrum. Later Tisza introduced the two-fluid 
model of superfluidity, further developped by Landau. At 96, he is still active at MIT. 

In the late twenties, von Neumann and Szilard even suggested that the university examina- 
tions could be replaced by such a contest [2J. 

The tradition still continues, and contributes to forming future generations of scientists 
PQ. The contests have been internationalized with the International Olympiads for High-School 
Students: first in Mathematics, and, since 1967, also in Physics. 



The Eotvos contest lasts 5 hours, and the use of all documents or tools is allowed. Three 
problems that require imagination and creative thinking rather than lexical knowledge are asked. 
They can, in some cases, lead to genuine research. A famous example is that of A. Csaszar (a 
distinguished topologist), who, as a young assistent, was called to survey the contest. While 
the high-school students were working, Csaszar figured out a generalization of the geometric 
problem given that year, and later published a paper on his findings. 

2 A problem of spacecraft landing 

Some of the physics problems deserve further thinking also. In 1969 - the year of the Moon 
landing, - for example, an exercise asked the following. A spacecraft moves on a circular trajec- 
tory of radius r = 4000 km around a planet of radius R = 3600 km. Then, for a short time, the 
rocket engins (directed oppositely to its motion), are switched on. This puts the spacecraft onto 
an elliptic trajectory, which touches the planet's surface at the opposite point of the trajectory. 
Which proportion of its kinetic energy went lost ? 

The problem can be solved using elementary methods alone 1 . It is, however, instructive to 
describe it also using different, more sophisticated, methods, which provide us with an insight 
into the problem of planetary motion. 

Prom a physical point of view, we have the following situation. As the engine works for a 
very short time, the position of the spacecraft does not change considerably. Owing to its loss of 
velocity, the gravitational attraction wins the race and pulls the circular trajectory more tight : 
the trajectory becomes an ellipse with major axis 2a = r + R. Our task amounts to comparing 
the kinetic energies of the two types of motions in the same (aphelium) point. 



Figure 1: Reducing the velocity puts the spacecraft onto an elliptic trajectory. The relative loss 
of energy is the numerical excentricity of the new orbit. 

The clue is to establish the following statement: The ratio of the [aphelium] kinetic energies 
of the two types of motion is the perihelium- distance divided by the major semi-axis of the 
elliptical motion, 
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x As the impact speed does not vanish, the problem is obviously not a realistic model of Moon landing - but 
this has not been its aim, neither. 



Then the answer to the originally asked question follows at once: 
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Now, as e = \{R + r) — R = |(r — R) is the excentricity, © is indeed e = e/a : the relative loss 
of energy is the numerical excentricity of the new orbit, i. e., the measure of the "flattening" of 
the circle. Let us observe that the answer only depends on the geometric dimensions. With the 
given numerical data, we find e = 1/19. 

Below we present several proofs of (^Q), ordered following their ever-increasing difficulty and 
required background knowledge. 

3 Demonstrations 



Proof I: using Kepler's laws alone . 

The first proof is elementary. According to the laws of circular motion, 

vJ = [M_ 

where / is Newton's constant. The [square of the] period is hence 
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This circular motion has kinetic energy 

fMm 
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Application of Kepler's third law to the elliptic motion yields the period of this latter, 

T 2 _{R + rf T _ 2tt /R + r\^ 2 (g) 
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The area of the ellipse is Trab, where b is the minor semi-axis, b 2 = V a 2 — e 2 = \J rR, the 
area is hence tt(R + r)V rR/2. The areal velocity, which is constant by Kepler's second law, is 
therefore 

nab I fMrR 
V = — = iWTr)- (7) 

At the aphelium v = \v a r, so that after slowing down, the velocity is 
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The corresponding kinetic energy is then 



E kin = ™i = f Mm R = E kin 2R (q) 
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which implies Q. 



Proof II: Using the conservation of the energy and of the angular momentum . 

Denoting the aphelium and the perihelium velocities by v a and v p , the conservation of the 
energy and of the angular momentum (divided by m) requires that 
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Eliminating the perihelium velocity yields once again the kinetic energy Q. 



(11) 



Proof III: Using the formula of the total energy of planetary motion . 

An important property of planetary motion E] is that the total energy only depends 
on the major axis, according to 

= JMrn 

2a V ; 

Then it follows from the energy conservation that, in any point of the trajectory, the velocity 
satisfies 

v 2 = fM ( I - I). (13) 
\r a J 

For the circular motion, a = r, and for the elliptic a = (r + R)/2, respectively. Plugging this 
into ()13j) . yields © and ©, respectively. 

Even more simply, observing that the change of the total energy is in fact that of the kinetic 
energy since the potential energy is unchanged, using (fT2*)) we have 

AE kin = AE tot = fM ( _±_ _L\= { fM^f_^R (14) 
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Writing here by © Hq™ in place of fM/2r yields © directly. 

It is worth noting that, at a point r of the trajectory, the ratio of the kinetic and the potential 
energies is, by E kin /EP ot = r j 2a - 1. We have therefore 

Ef n = -iE pot , resp. E* in = --^E pot ; (15) 

which yields again (pQ). 

Furthermore, while the total energy only depends on the major semi-axis, this is not so for 
the parts taken individually by the kinetic and the potential energies. According to (|15l) we 
have indeed 

E kin jjjtot jjjkin _J^°^ (16) 



Proof IV: From the radial equation 

The result can also be obtained from studying radial motion 2 . For an arbitrary central 
potential V, the problem can be reduced to one-dimensional motion with effective potential 
2 This was suggested to me by J. Kiirti. 



V e ff = V + 2mr' 2 wriere L i s the total angular momentum, cf. [3] p. 75. The radial velocity is 
therefore given by 
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where £7 = i? toi is the total energy. In the extremal points r a and r p the radial velocity vanishes, 
so that, for V = —fMm/r, 
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But the aphelium velocity is v a = L/mr a so that 



For a circular motion at r = r a , fM/r a = Vq, and therefore 
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consistently with our previous result. 

Proof V: Relation to Kepler's third law . 

Kepler's third law is related to the behaviour of the system with respect to scaling [3] : if 
some trajectory is dilated from the focus by A, and the time is dilated by A 3//2 , 

r^r' = Ar, t^t' = \ z ' 2 t, (20) 

yields again a possible trajectory. In those points which correspond to each other, both the 
kinetic and the potential energies [and hence also the total energy] are related as the inverse 
ratio of the geometric dimensions, 

f -A"'. (21) 

Let us now retract our original circular motion so that its radius equals to the major semi- 
axis of our elliptic motion above, i.e., consider the dilation by A = \{r + R)/r. By lf2~T|) the total 
energy [and consistently with (|12jl] is 
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This is, however, the same as the total energy of the elliptic motion, £q°' = E^, since the major 
semi-axes are equal. Hence once again 
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Then the result follows from (|16|). 

Let us stress that Kepler's third law did not suffice alone; we also needed the statement 
about the total energy. 



Proof VI: Using the Frenet formulas . 

It is worth describing the motion using the moving local frame introduced by Frenet [Jj. 
Then, for a trajectory of arbitrary shape, the normal component of the acceleration is v 2 /p 
where p is the radius of curvature i. e., the radius of the osculating circle [Zj. In an extremal 
point of the ellipse the accelaration is normal, and points toward the focus. Hence 



m — = Force 
P 

which generalizes the formula © of circular motion. For the circle p = r, so that 
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since the force is the same for both problems. We have hence proved: The ratio of the kinetic 
energies is identical to that of the radii of curvature. In the extremal points of the ellipse, 
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which implies again If!]). This confirms our intuition: decreasing the velocity increases the 
curvature. Using the explicit form, fMm/r 2 , of the force, (|22[) would allow us to calculate the 
velocity as 

' >f M , P__<v P 
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This is, however, not necessary for us: it was enough to know the geometric dimensions of the 
trajectory. 



Proof VII: Using the "Runge-Lenz" vector . 

A proof analogouos to that in II is obtained if we use the so called "Runge-Lenz" vector 

[311 El El 

K = mv x L — fMm r (25) 

where L = r x v is the conserved angular momentum; r denotes the unit vector carried by the 
radius vector drawn form the Earth's center to the spacecraft's position. 




Figure 2: The conserved Runge-Lenz vector is directed from the Earth' center towards the peri- 
helium position. 



Differentiating K-t with respect to time shows, using 

v (r • v) 1 1 
r = 5 — r = — (r x (v x r)) = — x CL/m) 

and the equation of motion, that K is a constant of the motion. The scalar product of K with 
L vanishes, so that K lies in the plane of the motion; it points from the focus to the perihelium 
point : K = Ke. Multiplying (|25|) with r yields the trajectory E] as 

r - E V =— e = — (26) 

l + ecosfl y fMm 2 fMm v ; 

where 6 is the angle between K and r. (j2fifl defines a conic with parameter p and numerical 
excentricity e. 

Returning to our initial problem, let us observe that in the extremal points 

K = mVpLe — fMme = —mv a Le + fMme, (27) 

where e is the unit vector directed from the center to the perihelium. The length of L is clearly 
L = mv p r p = mv a r a [cf. eliminating the perihelium velocity, 
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For circular motion r a = r p = r yielding for our elliptic motion r p = R, r a = r which 
provides us again with Q, the kinetic energy in the aphelium. Squaring (|27|) yields furthermore 

K 2 = f 2 M 2 m 2 + 2E tot L 2 . (29) 

Hence K = fMme which, together with (J27|) yields 

vl = l^(l-e). (30) 
Writing 2Eq iti for fMm/r provides us again with (0) or ©. 

Proof VIII. Using the hodograph 

Drawing the instantaneous velocity vector from a fixed point O of "velocity space" yields 
the hodograph. For planetary motion this is a circle [51 |Hj ■ The simplest proof of this statement 
is obtained if we multiply the angular momentum vectorially with the Runge-Lenz vector j^]. 
Developping the double vector product and using L • v = yields 

L x K = L 2 mv - fMmLu 

where u is the unit vector obtained by rotating r around the direction of L by 90° degrees in 
the counter-clockwise direction. As L x K = —LKj where j is the unit vector directed along 
the y axis of the coordinate plane. Writing here K = fMme the velocity vector is expressed as 

fM 



L v u - £ jj- ( 31 ) 

As the unit vector r turns around during the motion, so does also u (advanced with (90° 
degrees). The first term in (|31|) describes hence a circle of radius fM/L, whose center has been 
translated to C, situated on the y axis at distance —K/L = fMe/L below O. For a circular 
trajectory C = O and the hodograph becomes a circle around the origin with radius v 2 = fM/r. 
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Figure 3: The velocity vector of planetary motion describes a circle of radius fM/L and center at 
C, situated on the y axis at distance (OC = K/L = efM/L) in "velocity space". The aphelium 
point corresponds to the shortest distance from O, represented as the top of the circle. 

The velocity is the largest at the bottom of the circle, which corresponds to the perihelium 
point. The smallest velocity is obtained in turn in the aphelium, which is the top of the circle. 
Then u points vertically upwards, u a = j. Then the length of this smallest velocity is plainly 



which yields (|3U|) . Alternatively, we can write L = v a r a in (|31l) . 

Proof IX. Maxwell's generalization 

It is worth noting that our problem here can be readily generalized. J. C. Maxwell has in fact 
shown in a lovely little booklet that, in an arbitrary point r of the trajectory, the velocity is 



Maxwell' proof only uses Kepler's second law and some geometric properties of the ellipses. 
Let in fact the Sun be in one of the foci, S, let H be the other focus and draw perpendiculars 
SY and HZ from the foci to the tangent to the ellipse at a point P, see Fig. 4. 



smallest velocity is 
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(32) 




Figure 4: Maxwell's geometric derivation. 



By Kepler's second law, the areal velocity is constant and given by twice the area of the 
ellipse, nab, divided by the revolution time, T, 

w 2nab 2nab 1 

SY • v = — =>• v = — , (33) 

T T sr 

On the other hand, a known geometric property of the ellipse says that the half minor axis, 



b, is the mean of SY and HZ, 



SY • HZ = b . 



Finally, the "optical property" [known to Kepler] says that all light rays emitted from one 
focus are reflected by ellipse to the other focus. It follows that the angles HPZ and SPY are 
equal; the triangles SPY and HPS are therefore similar. Hence 



SP : PH = SY : ZH. 
Multiplication of the two relations yields 

W 2 = b 2 SP 



PH 2a -r 
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since SP + PH = 2a, the major axis of an ellipse. Inserting SY into the square of (|33|) yields 
Maxwell's formula (f52|) . 

For a circular orbit in particular, v 2 = An 2 r 2 /Tq. (In general, i. e. with the exception of the 
extremal points, the velocities vo and v have different directions.) Hence 

v 2 _ a 2 T 2 2a - r 
Vq r 2 T 2 r 

Using here Kepler's III law, we end up with the following simple expression: 

v 2 2a -r 
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In words: At any point of the orbit, the square of the velocity to that of the circular motion 
through the same point is as the distance from the other focus to the semi major axis ! For 
r = r a we recover in particular our previous result 

Proof X: From Newton's Principia 

Last but not least, let us point out that the problem discussed here has actually been solved 
by Newton over 300 years ago ! Proposition 16. Theorem 8. Corollary 3. of his "Principia" [5] 
says in fact that "the velocity in a conic, at the greatest or least distance from the focus, is to the 
velocity with which the body would move in a circle, at the same distance from the center, as the 
square root of the principal latus rectum is to the square root of twice that distance." Translated 
into our notations, 

v 2 a :v 2 =l:2r. (35) 

But the half the "principal latus rectum" is the harmonic mean of the periapsis and the apoapsis 
of the ellipse ^U] (and also the parameter p = b 2 /a), 
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Inserting into (|35jl the result (^Q) is recovered once again. 

Newton's original proof comes by a series of elementary geometric demonstrations using 
hardly more then the properties of similar triangles 9: ; . The reader is referred to Newton's 
Principia for details. Let us remark, though, that (|34j) provides a quick proof to Newton's 
statement (|35|). since the prodcuct of r = r a = a(l + e) and r p = a(l — e) is r a r p = a 2 (l — e 2 ) = 
al/2. 

4 Discussion 

Let us summarize our various approaches. Our first proof only used Kepler's laws specific for the 
planetary motion, and suits perfectly to a high-school student. The second and seventh proof 
is based on conservation laws; the second uses that of the energy and the angular momentum, 
and the last the Runge-Lenz vector. This is early 19th century physics: the vector (|25j) was in 
fact introduced by Laplace in 1799, in his Traite de Mecanique Celeste [HHUE]- 
Proof IV is based on the radial equation. 

II, using high-school knowledge only, would clearly work for any conservative force, while 
VII is related to the "hidden" symmetry of of the Kepler problem. Although some of the results 
could be obtained by freshmen, this approach is not in general taught at the university. It leads 
to a group theoretical treatment of planetary motion !6i. For example, L • K = and (|29|) are 
the classical counterparts of the Casimir relations of the SO (4) dynamical symmetry used by 
Pauli to determine the spectrum of the hydrogen atom [H] . 

Ill and V are based on the formula (fl"2|) of the total energy, discussed by university textbooks 
[HI V is linked to the scaling property which yields in fact Kepler's third law (5]. 

Proof V uses the general framework of co-moving coordinates called the Frenet formula? 
[Zj (late 19th century), which makes part of regular university courses on mechanics and/or 
differential geometry. It can be applied to any central force problem: the reader is invited to 
work out what happens, e, g., for a harmonic force F = —kr (when the trajectories are again 
ellipses.) 

Proof VIII is based on the hodograph, also used by Feynman in his geometric approach to 
planetary motion jSj. 

The problem is readily generalized following Maxwell (Proof IX) , along the lines of Proof 
I and adding some knowledge of geometry. Remarkably, this general answer might have been 
found by Kepler in the early 17th century ! 

Our final proof X uses a statement made by Newton in his Principia, also obtained by 
elementary geometry. 

Interestingly, the very first question of the very first Eotvos physics contest (that of Szilard 
in 1916), was on planetary motion: Let us consider a miniaturized model of the solar system, 
where all distances are scaled as 1 : 15 • 10 10 . How long would last a year ? The answer is that it 
would remain a year, since, assuming that all densities remain constant, the Sun's mass would 
also decrease as the cube of its dimension. But this is precisely the condition, as seen from the 
scaling argument in Proof V • 
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